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ABSTRACT ■ 

An approximate boundary for a finite sample sequential decision 
process is presented (without derivation). By means of Computer Simu- 
lation, properties of the process are checked; the power of the test is 
thus determined for selected alternative hypotheses. 
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I. INTRODUCTION 



The main purpose of this paper is to study the properties of a 
mathematically derived sequential sampling procedure. The derivation 
involves the assumption that a certain parameter N, identified as a 
sanple size, was "large." Our purpose was to make computer calculations 
that would reveal how the sampling procedure might behave for finite, 
or even rather small, N. Other calculations would also indicate the 
power characteristics of the test, that is, the probability of rejecting 
a null hypothesis. The null hypothesis, in this case, was that a partic- 
ular probability distribution describes the sample . By experimental 
sampling methods we attempted to estimate the power that the procedure 
possessed to detect departures of various sorts from the null hypothesis. 

The stochastic process used to obtain the sequential procedure might 
also be used for other purposes in operations research. For example, 
the sequential procedure might be used to describe arrivals in transitory 
queueing situations. Some of these situations are briefly described at 
the end of this paper. 
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II. A SEQUENTIAL SAMPLING PROBLEM 



Suppose, for example, that a supplier of transistors asserted that 
each of his transistors had an in-service lifetime T, where T was a 
random variable distributed exponentially; 

r 1 - e , t ^ 0 

F(t) 

L 0 , t < 0 . 

In order to test that assertion or hypothesis, a sample of N tran- 
sistors was selected and placed in service simultaneously at time t = 0. 
On the average, the number of transistors that failed by time t >0 was 
NF(t), i.e.- if denoted the number of transistors that failed by 

time t, then 

E [A^j(tj] = uF(t) 
and 

Varfjj(t)] = NF(t) [l-F(t)] , 

since the transistors i,;ere presumed to fail, or die, independently. 

If one waited until all N transistors had failed, the rate coiild be 
estimated and compared to A hypothesis test can be conducted by 

waiting until r out of N transistors failed (r<N); this was shown by 
Epstein and Sobel [13 . 

The present approach to this problem was to establish a boiindary^ 
Bjj(t), with the property that if >= i.e. if the supplier was en- 

tirely honest, then Ajj(t)< Bjg(t) for all t, with a known, specifiable, 
probability denoted by 1-6. For instance, suppose 6 = 0.05« Then, if 
^ the probability that the accumulated numbers of failures by 
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time t vrould ever exceed the boundary is five per cent. Hence, if 
remains below the boundary, £ A with ninety-five per cent confidence. 

The detailed derivation of the boundary is given by Gaver, Lehoczky, 
and Perlas in an unpublished technical report that is in preparation. The 
result is quoted below. 

The Boundary 

If F(t) is sufficiently smooth (possesses a continuous density) and 
if 

Bjj(t) ,= NF(t) +VN {c<F(t) [l-F(t)]] , (2.1) 

then the probability that •^('t) < ^ is approximately 

-2o(p 

1 -e j the approximation irproving as N becomes large (N—»CD). Here, 
c< and p are arbitrary positive numbers . 

A graph of a possible failure arrival pattern is a step function with 
unit increases at time points that are the order statistics of F(t): 

T(i ), T^ 2 )^ ••• where ) is the smallest of the N failure times, T^ 2 ) 
is the next smallest, ... T^j^j is the kth smallest. Then, a particular 



sample nay be represented as 


fo].lows 


2 


< 

c+* 

(1 

O 


for 




= 1 


for 


' ’■(2) 


= 2 


for 

• • • 




= k 


for 

• • • 


’■(k) ' * ‘\k+1) 


= N 


for 


’’(K) ^ * ’ 


and graphically, we have the 


following picture. 



( 2 . 2 ) 
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In Figure 1, a crossing of the boundary occurs at T^. We would conclude 

that the null hypothesis is not satisfied if such an event occurs. Notice 

that the parameters ot and 'p may be selected at will. It is interesting 

_2«<'R -2k 

to arrange that o< p = k, a constant such that e •^ = e = 0.0^, 
for example^ and then to attempt to select o< ( or 'p ) so that the 
boundarj’’ will be crossed soon if a non-null distribution of specific 
sort is to be detected [2]. This problem has not been extensively in- 
vestigatedj but the expected time to hit the boundary has been tabulated 
for certain of the problems studied. 
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III. USE OF THE COMPUTER SIMULATION 



Several facts about the boundary %(t) presented in equation (2.1 ) 
are worthy of notice . 

-2«B 

a. The prescribed probability of crossing the boundary, e ' , 
is not exact for finite N. That is, use of this probability to assess 
the probability of crossing is an approximation which is, in fact, quite 
analogous to the use of the normal distribution as an approximation to 
the distribution of a sum of independent random variables. The mathema- 
tical derivation does not indicate whether the actual probability of 
crossing is, perhaps, consistently lower or higher than the actual true 
probability. More complex mathematical methods may provide useful in- 
formation of this type, but a straightforward approach is that of simu- 
lation or synthetic sampling. By this iiBans, it can be seen how increases 

/ 

in N improve the quality of our approximation, as is suggested by the 
mathemtical derivation . 

b. The relationship of the true crossing probability to the 

approximate crossing probability, for given N, may possibly be 

influenced by the choice of o< and . Setting p = k, where k is a 
constant delivering a desired probability of crossing (for k~1 .ii9), 
o< (and hence p ) may be adjusted so as to obtain good agreement between 
the theory (e = e > ) and the true probability for various values 

of N — in particular for small N. In this paper, an attempt was made 

to study this problem by simulation. 

In order to study the properties of the boundary Bjj(t), the following 
simulation was carried out. 
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A. SIMULATION LOGIC 



1. The null hypothesis distribution was chosen and called FQ(t). 
The formula used to calculate the boundary values was as follows: 

Bj^(t) = NFQ(t) +'TN • (3-1) 

An exponential distribution is the null hypothesis, as in the illustra- 
tive problem; hence, for the initial calculations, the equation 

%(t) = N(l-e'^) +fN (3-2) 

was used. 

2. A sairpling distribution was selected. Initially, it was the 
null hypothesis distribution, so that the accuracy of the theoretical 
estimate of the probability of crossing the boundary could be checked. 
Later, other distributions were tried in order to determine the power 
for departures from the exponential distribution. 

3. A sample of size N was taken and arranged in order: 






(3-3) 



This ordered sample was conposed of the actual numbers obtained as 
realizations of the random variables. 








(3.U) 



The details on how the random samples were conducted will be presented 
shortly . 

1^. At each failure time, say starting with k = 1, then 2, 3, 

..., = k was compared with if For k = 

1, 2, ..., the process continued to k + 1 and was checked again. If at 
k--, the event k* ^ occurred for the first time,' a boundary 



crossing was recorded at the time t^ 5 the test was stopped and 

another simulation begun. 

5 . Step four^ above j was independently repeated NR times (NR means 
the niimber of replications). The probability of crossing f^^ was esti- 
mated by the fraction 



nvimber of crossings. 
% " NR 



(3.5) 



Clearly, NRfjj was a realization of a binoraially distributed random 
variable. An approximate idea of the variability of the estimate of 
the probability of crossing is given by 



NR 



NR 



(3.6) 



An approximate confidence limit for the true probability of crossing 

^ 2 

came from assuming that f.. is approximately normal with variance 

NR 

B. SIMULATION PROGRAM 

The simulation procedure just described was written in FORTRAN com- 
puter programming language and r^un on the IBI-I 300, Model 67 Computer at 
the Naval Postgraduate School. The program was written to facilitate 
changes of distributions, boundary parameters ( « and "p ), and numbers 
of replications (NR). It was also constructed so that the saii^jling 
process could be stopped after a fixed time period, i.e.- before all of 
the order statistics had been recorded. 

The pseudo random number generator employed to obtain realizations 
of uniform random numbers was a library subroutine called RANDOM [3 and U] ; 
it is based on a Lehmer multiplicative scheme. The uniform realizations 
were then transformed into other random numbers by a probability integral 
transformation . 
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A program listing and sample output has been included at the end 
of this paper. In order to generate the sequence of order statistics 
(interpreted as failure or arrival times), the random sample of size N 
was first created. Then, a SORT routine was employed to put the times 
in order. A more efficient procedure would have been to generate the 
ordered increments sequentially; for example, the time between and 



(k+1 ) 



in exponential sampling is also exponential with a mean propor- 
1 



tional to One program has been written in this vein, but production 

runs for the present investigation have largely been carried out by use 
of the SORT routine 
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IV. SIMULATION EXPERIMENTS 



Several siitiulation experiments have been carried out to develop an 
\inders landing of the boundary crossing process. 

A. EXPERIMENT 1 

This experiment explores the crossing probability under the null 
hypothesis ( = ^q) and under selected alternatives, for various sample 

sizes (N)j and for different boundary parameters o( and"^ , but with 
<< = "P . The computed probability of a crossing fj^ was calculated by 
the fraction ( 3 * 5 ) mentioned previously. The average number of arrivals 
up to a crossing (JT(t)) was calculated by summing the number of arrivals 
(Ajj(t)) that had occurred up to the time of crossing, over the two thou- 
sand replications and dividing by the number of crossings (NC) recorded. 
The average crossing time (C(t)) was calculated by adding the arrival 
times, which generated crossings, then dividing that sum by the number 
of crossings. The results appear in Table I (theoretical crossing prob- 
ability of five per cent), and Table II (theoretical crossing probability 
of ten per cent). These Tables are based on two thousand replications 
(NR = 2000). 

When the boundary was not crossed for a sample size N, all values for 
that N appear as 0. In such a situation, the zero (O) is intended as a 
meaningless place holder rather than the unique number. 
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TABLE I. 



Probability of a Crossing 







«<= ■p = 1 .22388 










EXPON ( >» = .5) 






N 


NC 


A(t) 


C(t) 


£n 


U 


3 


3.33 


.133 


0.15 


6 


2 


5.00 


.300 


0.10 


8 


7 


6.U2 


• U28 


0.35 


10 


1 


8.00 


.523 


0.05 


25 


0 


0.00 


.000 


0.00 


50 


0 


0.00 


.000 


0.00 


75 


0 


0.00 


.000 


0.00 



N 


NC 


EXPON ( A = 1 ) 
A(t) 


c(t) 


!n 


h 


U6 


3.63 


.287 


2.30 


6 


69 


5-07 


.391 


3-U5 


8 


66 


6.30 


.U02 


3.30 


10 


70 


7.37 


.Ul6 


3.50 


25 


73 


15.19 


.U61 


3-65 


50 


85 


27. 5U 


•5oU 


U.25 


75 


lOli 


U2.l43 


.586 


5.20 







EXPON ( A = 1 .25) 






N 


NC 


A(t) 


C(t) 


%■ 


h 


100 


3.66 


.298 


5.00 


6 


liil 


5.11 


.381 


7.05 


8 


165 


6.33 


•U111 


8.25 


10 


197 


7.U8 


.U28 


9.8U 


25 


372 


15.95 


.513 


18.59 


50 


61;U 


28. 6U 


•5U2 


32.19 


75 


839 


U1 .09 


.560 


.95 
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TABLE I. (continued) 



EXPON ( ^ = 1 . 5 ) 



N 


NC 


A(t) 


c(t) 


_ 5 l 


h 


171 ; 


3.68 


.303 


8.70 


6 


253 


5.03 


.362 


12 . 61 ; 


8 


331 


6.26 


.395 


16 . 5 U 


10 


liOi; 


7.i;8 


.Ii 26 


20.19 


25 


896 


15.67 


.1;98 


hh.79 


50 


11;61 


25. 5 U 


.li 70 


73. 01 ; 


75 


1770 


35.59 


. 1;32 


88 . 1;9 







EXPON ( A = 1 . 75 ) 






N 


NC 


A(t) 


C(t) 




1; 


262 


3.70 


.301 


13.10 


6 


389 


1;.97 


. 3 l ;9 


19 . 1;5 


8 


525 


6.19 


.381 


26.21; 


10 


672 


7 . 1;3 


.1;18 


33.59 


25 


1391 


11 ;. 80 


. 1;37 


69 . 51 ; 


50 


1893 


23.31 


.362 


9 l;. 61 ; 


75 


1985 


28.30 


. 271 ; 


99. 21 ; 







EXPON ( X = 2 ) 






N 


NC 


A(t) 


c(t) 




1 ; 


360 


3.69 


.296 


18 . 1;5 


6 


552 


U .98 


. 3 l ;6 


27.59 


8 


7 U 1 


6.15 


.372 


37. oU 


10 


929 


7.31 


.397 


1 ; 6 . 1 ; 1 ; 


25 


1736 


13.95 


.382 


86.79 


50 


1985 


19.86 


.256 


99. 21 ; 


75 


2000 


23.65 


.190 


100.00 
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TABLE II. 



^0% Probability of a Crossing 
^ .07298 

EXPON ( >> = .^) 



N 


NC 


A(t) 


C(t) 




h 


15 


3.33 


.300 


0.75 


6 


9 


5.00 


.1^67 


0.1*5 


8 


11 


5-72 


.31*1 


0.55 


10 


5 


6.80 


.liOO 


0.25 


25 


0 


0.00 


.000 


0.00 


50 


0 


0.00 


.000 


0.00 


75 


0 


0.00 


.000 


0.00 



EXPON ( > = 1 ) 



N 


NC 


A(t) 


c(t) 




1* 


127 


3.1*9 


.388 


6.35 


6 


151 


1*.90 


.1*28 


7.51* 


8 


151* 


5.92 


.1*27 


7.70 


10 


156 


7.03 


.1*1*3 


7.79 


25 


160 


1ii.86 


.1*95 


8.00 


50 


177 


26.71* 


.522 


8.85 


75 


190 


39.06 


.551 


9.50 



EXPON ( > = 1 .25) 



N 


NC 


A(t) 


C(t) 


5l 


1* 


21*1 


3.51 


.31*5 


12.05 


6 


289 


1*.72 


.381 


11*.1*5 


8 


331* 


5.98 


.1*25 • 


16.70 


10 


371* 


7.10 


.1*1*6 


18.70 


25 


627 


15.11* 


.517 


31 .31* 


50 


91*5 


26.86 


.527 


1*7. 2li 


75 


1168 


38.66 


.51*3 


58.39 






TABLE II. (continued) 
EXPON ( A = 1 ,g) 



N 


NC 


A(t) 


c(t) 


£n 


h 


360 


3.51i 


.31|5 


17.99 


6 


1;60 


1;.72 


.376 


22.79 


8 


571 


5.85 


.1;05 


28.51; 


10 


67 k 


7.0li 


.1;39 


33.69 


25 


1197 


U.li5 


. U 72 


59.81; 


50 


1682 


21;. 00 


• h 2 h 


81;. 09 


75 


1907 


31.37 


.372 


95.3I; 







EXPON ( > = 1 .75) 






N 


NC 


A(t) 


C(t) 




h 


U 9 S 


3.53 


.337 


21;. 89 


6 


660 


h .7 h 


.381 


32.99 


8 


823 


5.81; 


.399 


1;1.11; 


10 


972 


6.88 


.lilO 


1;8. 59 


25 


1657 


13.59 




82.81; 


50 


1957 


20.36 


.307 


97 . Sh 


75 


1998 


2l;.l;6 


.230 


99.89 







EXPON ( > = 2 ) 






N 


NC 


A(t) 


C(t) 


5 l 


1; 


617 


3.51 


.326 


30.81; 


6 


868 


1;.71 


.370 


1+3.39 


8 


1083 


5.77 


.387 


51;.ll; 


10 


I2l;7 


6.71; 


.386 


62.3 k 


25 


1878 


12.1;0 


.337 


93.89 


50 


1997 


17.22 


.211; 


99.81; 


75 


2000 


20.51 


.160 


100.00 
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B. DISCUSSION 



The actixal probability of crossing appeared to be consistently smaller 
than the theoretical probability under the null hypothesis (A = 1 ). For 
N of fifty, and over, the agreement appeared to be reasonably good. More- 
over, the probability of rejecting the null hypothesis (power) was quite 
high when A = 1 and N was at least twenty -five. This would not be an 
unreasonable sample size for the transistor testing illustration, but 
would undoubtedly by large if the procedure were applied to check entire 
systems . 

Notice that an empirically generated correction could be derived to 
adjust the size of the null hypothesis probability. Consider, for ex- 
ample, the N = li, A = 1 .2^ element of the nominal five per cent table: 
the crossing per cent (estimated true probability of crossing) is nearly 
five per cent. Under the null hypothesis, the simulation crossing prob- 
ability was computed at 2.3 per cent. Noxj, if the null hypothesis pre- 

-1 

vailed, and if each failure time were multiplied by (1 .25) = 0.8, then 

five percent boundary crossings vfould be realized. On the other hand, 
if A > 1 , and each failure tiine was multiplied by 0.8, but still applied 
against the theoretical derived boundary for A =‘ 1 , the power of the 
procedure would be further enhanced. This empirical adjustment has been 
left for more extensive development, e.g. — for different sample sizes — 
in the future. It seems clear that an alternate approach involves an 
adjustment of the boundary, itself. 

C. EXPERIMENT 2. 

Using the same procedure as in Experiment 1 , we explore the effect 
/ 

of o< ^ p , but with ■p = constant. We tried the follovring four com- 
binations; (1 ) = 1 , (2) ^ = 1 .5'=< , (3) = 2f , (I;) ^ = 2o^ . 

The results are summarized in the next tables. 
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TABLE III 



Probability of a Crossing 
= 1 .1i9893 f = .99929 



EXFOM ( >> =.^) 



N 


NC 


A(t) 


c(t) 


f 

Ji 


h 


2 


i^.OO 


.156 


.10 


6 


1 


6.00 


.1|76 


.05 


8 


3 


7.33 


•561i 


.15 


10 


0 


0.00 


.000 


.00 


25 


0 


0.00 


.000 


.00 


50 


0 


0.00 


.000 


.00 


75 


0 


0.00 


.000 


.00 



EXFON ( ^ = 1 ) 



N 


NC 


A(t) 


c(t) 




1| 


17 


1^.00 


.297 


.85 


6 


36 


5.6i 


.ii6l 


1 .80 


8 


hh 


6.95 


.li92 


2.20 


10 


h9 


8.20 


.5oh 


2.U5 


25 


h9 


17.28 


.589 


2.1;5 


50 


70 


32.91 


.697 


3.50 


75 


91 


h7-h6 


.706 


1^.55 



EXPON ( ^= 1 .25) 



N 


NC 


• A(t) 


c(t) 


In 


h 


36 


U.oo 


.309 


1 .80 


6 


97 


5.59 


.ii57 


ii.85 


8 


118 


7.00 


.li93 


5.90 


10 


1ii7 


B.3h 


.529 


7.35 


25 


3liO 


18.20 


.668 


16.99 


50 


636 


32.99 


.697 


31 .79 


75 


871 


h6.$h 


.692 


li3.5ii 
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TABLE III (continued) 
EKPON ( = 1 .g) 



N 


NO 


A(t) 


c(t) 


i 


\x 


76 


1;.00 


.316 


3.80 


6 


161 


5.59 


.1;28 


8.01; 


8 


232 


7.03 


.1;95 


11 .60 


10 


321; 


8.36 


.531 


16.19 


25 


857 


17.i;7 


.603 


1;2.81; 


50 


11;72 


29.77 


.562 


73.59 


75 


1799 


l;0.3i; 


.521 


89.91; 



EXPON ( >> = 1 . 1 ^) 



N 


NO 


A(t) 


c(t) 




1; 


126 


U.oo 


.311; 


6.30 


6 


25I; 


5.55 


.i;l5 


12.70 


8 


1;01 


6.91 


.1;70 


20.01; 


10 


560 


8.27 


.511 


27.99 


25 


1381 


16.67 


.538 


69.01; 


50 


1897 


26.1;2 


.1;36 


91;.81; 


75 


1990 


32.57 


.331; 


99. h9 






K(PON ( >> = 2) 







N 


NO 


A(t) 


c(t) 


!n 


1; 


189 


1;.00 


.310 


9.1U; 


6 


1;02 


5.50 


.1;09 


20.09 


8 


591; 


6.85 


.1;51; 


29.70 


10 


787 


8.11; 


.1;81 


39. 31; 


25 


1733 


15.70 


.1;62 


86.61; 


50 


1987 


22.91 


.312 


99. 31; 


75 


2000 


27.50 


.229 


100.00 
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TABLE IV 



Probability of a Crossing 
^ = .99929 f = 1 M 93 



EXPON ( ^ = .5) 



N 


NC 


A(t) 


G(t) 


5l 


li 


10 


3.10 


.203 


0.50 


6 


h 


It. 50 


.273 


0.20 


8 


10 


5.70 


.331 


0.50 


10 


3 


5.66 


.2lt0 


0.15 


25 


0 


0.00 


.000 


0.00 


50 


0 


0.00 


.000 


0.00 


75 


0 


0.00 


.000 


0.00 



EXPON ( >> = 1 ) 



N 


NC 


A(t) 


C(t) 




It 


80 


3.30 


.2li1 


il.OO 


6 


99 




.301 


U.95 


8 


90 


5. lilt 


.291 


U.50 


10 


101 


6.27 


.305 


5.05 


25 


91 


12.68 


.331 


il.55 


50 


99 


23.08 


• 37ii 


li.95 


75 


117 


3li.13 


.Itll 


5.85 






EXPON ( A = 1 .25) 







N 


NC 


A(t) 


c(t) 


fa 


h 


162 


3.38 


.263 


8.10 


6 


198 


it. 52 


.301 


9.90 


8 


21 It 


5.61 


.323 


10.70 


10 


2li2 


6.57 


.33ii 


12.10 


25 


392 


13.78 


.392 


19.59 


50 


628 


2lt.36 


.li13 


31 .39 


75 


797 


35.16 




39.8it 
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TABLE 17 (continued) 










EXPON ( = 1 .5) 






N 


NO 


A(t) 


c(t) 


Si 




2li8 


3.U2 


.268 


12.39 


6 


33li 


li.li? 


.291 


16.70 


8 


liOO 


5.59 


.319 


19.99 


10 


U55 


6.58 


.333 


22.71 


25 


889 


13.77 


.396 


hh-hS 


50 


lli06 


23. UU 


.386 


70.29 


75 


1720 


3l.2li 


.357 


85.99 



EXPON ( >i = ^ .7g) 



N 


NO 


A(t) 


C(t) 


Si 


h 


3 h 6 


3.ii3 


.261i 


17.29 


6 


hn 


l;.ii7 


.287 


23.8U 


8 


588 


5.55 


.313 


29.39 


10 


729 


6.60 


.335 


36. Ii5 


25 


1368 


13.16 


.361 


68.39 


50 


1867 


20.58 


.303 


93. 3U 


75 


1979 


2li.88 


.235 


98.95 







EXPON ( X= 2) 






N 


NO 


A(t) 


c(t) 




k 


U66 


3.ii5 


.268 


23.29 


6 


6 hh 


li.51 


.293 


32.19 


8 


806 


5.53 


.307 


I1O.29 


10 


977 


6.53 


.323 


li8.81i 


25 


17U 


12.37 


.316 


85.70 


50 


1980 


17.39 


.216 


98.99 


75 


not run 
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TABLE V 



Probability of a Crossing 

o<= 1.73082 p= . 865 U 1 







EXPON ( > = .3) 






N 


NC 


A(t) 


c(t) 




1; 


2 


I4.OO 


.136 


0.10 


6 


1 


6.00 


.h76 


0.03 


8 


1 


8.00 


.6h9 


0.03 


10 


0 


0.00 


.000 


0.00 


23 


0 


0.00 


.000 


0.00 


50 


0 


0.00 


.000 


0.00 


73 


0 


0.00 


.000 


0.00 







EXPON ( ^ = 1 ) 






N 


NC 


A(t) 


C(t) 


5l 


h 


h 


h.oo 


.169 


0.20 


6 


17 


3.88 




0.83 


8 


28 


7.30 


.36ij 


1 .iiO 


10 


3k 


8.76 


.367 


1 .70 


23 


39 


18.92 


.703 


1 .93 


30 


39 


36.86 


.838 


2.93 


73 


83 


32.07 


• 81;3 


U.23 







EXPON ( = 1 .23) 






N 


NC 


A(t) 


C(t) 


fk 


h 


11 


k.OO 


.207 


0.33 


6 


k8 


3.87 


.ii6o 


2.i;0 


8 


70 


7 . ho 


.338 


3.30 


10 


103 


8.8h 


.383 


3.13 


23 


308 


19.811 


.803 


I3.ii0 


30 


612 


33.67 


.808 


30.39 


73 


831 


30.30 


.803 


U2.3U 
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f 




TABLE V (continued) 
EXPON ( > = 1 .^) 



N 


NO 


A(t) 


c(t) 




li 


17 


1;.00 


.198 


0.85 


6 


91; 


5.86 


.1;1;1 


1;.70 


8 


161 


7.i;0 


.531; 


8. Oh 


10 


253 


8.89 


.595 


I2.6h 


25 


801; 


18.82 


.691 


ho. 19 


50 


1 1;63 


32.1;2 


.651 


73. 1h 


75 


1802 


1;3.71 


.589 


90.09 



EXPON ( > = 1 .7^) 



N 


NO 


A(t) 


c(t) 


% 


h 


28 


h.oo 


.20h 


1 .ho 


6 


169 


5.81 


.h27 


8.h5 


8 


286 


7.33 


.507 


lh.29 


10 


hl;7 


8.83 


.571 


22. 3h 


25 


1327 


18.03 


.620 


66. 3h 


50 


1895 


28.73 


.h92 


9h.7h 


75 


1987 


35.92 


.381 


99. 3h 







EXPON ( X = 2) 






N 


NC 


A(t) 


C(t) 




h 


5h 


h.oo 


.216 


2.70 


6 


239 


5. 81 


.hi 8 


11.95 


8 


hh7 


7.29 


.h97 


22.3h 


10 


666 


8.70 


.538 


33.29 


25 


I72h 


17.12 


.535 


86.20 


50 


1985 


25.28 


.358 


99. 2h 


75 


2000 


30. 7h 


.265 


100.00 
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TABLE VI 



Probability of a Crossing 
o(= .86^1i1 p = 1 .73082 



EXPON ( = .g) 



N 


NC 


A(t) 


C(t) 




k 


III 


2.78 


.151 


0.70 


6 


6 


ii.l6 


.258 


0.30 


8 


. 11 


5.5I1 


• 3U 


0.55 


10 


5 


1|.80 


.176 


0.25 


25 


0 


0.00 


.000 


0.00 


50 


1 


13.00 


.130 


0.05 


75 


1 


18.00 


.130 


0.05 



EXPON ( A = 1 ) 



N 


NC 


A(t) 


C(t) 


% 


li 


91 


3.02 


.203 


ii.85 


6 


117 


3.95 


.230 


5.85 


8 


106 


5.00 


.25U 


5.30 


10 


115 


5.60 


.239 


5.75 


25 


105 


11 .Ii8 


.281 


5.25 


50 


119 


20.31 


.301* 


5.95 


75 


122 


30.05 


.336 


6.10 



EXPON ( > = 1 .2$) 



N 


NC 


■ A(t) 


C(t) 






192 


3.10 


.225 


9.60 


6 


228 


I|.13 


.253 


11 .I|0 


8 


2li5 


5.10 


.27I1 


I2.2l| 


10 


265 


5.98 


.279 


13.25 


25 


1|00 


12.1|2 


.328 


19.99 


50 


600 


21 .I|2 


.337 


29.99 


75 


777 


31 .81 


.373 


38.8I1 
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TABLE VI (continued) 
EXPON ( >>=1-5) 



N 


NC 


A(t) 


C(t) 




ii 


299 


3.15 


.232 


li;.95 


6 


361 


it. 11; 


.253 


18.01; 


8 


ii22 


5.08 


.269 


21 .09 


10 


ii8l 


6.08 


.287 


21;. Oi; 


25 


859 


12.28 


.323 


lt2. 95 


50 


1363 


21 .33 


.331 


68.11; 


75 


1679 


28.53 


.312 


83. 91; 







EXPON ( > = 1 .75) 






N 


NO 


A(t) 


C(t) 


Si 


h 


399 


3.15 


.226 


19. 91; 


6 


501; 


1;.17 


.251; 


25.20 


8 


613 


5.07 


.263 


30.61; 


10 


728 


6.03 


.281 


36.39 


25 


1320 


11.96 


.309 


65.99 


50 


181;3 


19-00 


.271 


92.11; 


75 


1972 


22.88 


.212 


98.59 







e:':fon ( >> = 2) 






N 


NC 


A(t) 


c(t) 




it 


526 


3.15 


.229 


26.29 


6 


682 


1;.15 


.250 


3i;.09 


8 


837 


5.06 


.262 


i;1 .81; 


10 


991 


6.01 


.280 


1;9.51; 


25 


1703 


11 .36 


.278 


85.11; 


50 


1976 


15.98 


.191; 


98.79 


75 


not run 
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D. DISCUSSION 

The tendency of the theory to understate the crossing probability, 
noted ^ is even more pronounced when o(= 1 • However, when 

■p = 1 (TABLE IV) the tinders tat ement error nearly vanishes, and the 

test is very close to being exact (i.e.- crossings occur -wath probability 
about five percent) for all values of N checked. According to TABLE VI, 
the probability of rejection is slightly overstated when p = 2*^ , for 
all but N = ii. To within limits of sampling accuracy, it appears that 
p equalling 1 .6$o( will give a veiy nearly correct answer at 6 = 0.05. 
Further sampling will pin this down more precisely. A theoretical ex- 
planation V70uld be most interesting, for the asymptotic (large N) mathe- 
matics detects no difference in situations for which o(,p = constant. 

The next set of experiments vxas conducted in order to investigate 
our procedure's ability to detect departures from the null (exponential, 

A = 1 ) hypothesis. 

E. EXPERIIENT 3. 

Exploration of the crossing probability when the boundary reflected 
the null hypothesis (exponential vri.th ^ = )»q = 1, , nominal 

crossing probability five percent), but samples were dra-wn from a special 
long-tailed distribution with a Pareto-like appearance : 

^ , ax 

^p^^^ " 1+ax " ^ ^ 

In order to give a fair comparison, the parameter a was chosen so as to 
match the medians of F and the null distribution; thus a = (log„2) \ 

A rationale for the choice of (i;.1 ) is that such a distribution 
appears very similar to the exponential, but has a decreasing failure 
rate and a bias tovjard early failures. The results are recorded in 
TABLE VII. 
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TABLE VII 



5^ Probability of a Crossing 

o<= p = 1 .22388 





* 


PARETO 






N 


NC 


A(t) 


c(t) 


% 


h 


79 


3.50 


.253 


3.95 


6 


97 


I1.82 


.319 


U.85 


8 


91 


5. 61 


.271I 


i^.55 


10 


lOli 


6.70 


.309 


5.20 


25 


133 


12.30 


.277 


6 . 6 h 


50 


209 


20.00 


.257 


10. Ii5 


75 


218 


26.65 


.2l|2 


10.90 


F. 


DISCUSSION. 









Our sampling experiirent indicates that the sequential procedure imder 
study has very little power against the "Pareto tail" alternative (U.1 ) 
under the particular choice of parameters, at least for small N. Actually, 
the probability of rejection is somewhat larger than for the correspond- 
ing exponential. It seems likely that the relative prevalence of small 
obser^/ations accounts for the increased crossing probability, and also 
for the relatively short expected time to cross, given that a crossing 
occurred . 



G. EXPERIMENT U. 

The procedure used in Experiment 3 was repeated, except that samples 
were dravm from a shorter-tailed (Weibull, with increasing hazard rate) 
distribution: 



F„(x) 




(I1.2) 



Again, the parameter a was fitted by matching medians of the exponential 
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-1 

to by coincidence, a = (logg2) once again — an accident of 

parameterization. The results are recorded in TABLE VIII. 



TABLE VIII 

Probability of a Crossing 
e^= P = 1 .22388 



N 


NC 


A(t) 


c(t) 




h 


13 


ii.OO 


• ii 25 


0.65 


6 


38 


5.89 


.610 


1 .90 


8 


58 


7.72 


.702 


2.90 


10 


91 


9 .ii 6 


.778 


ii.55 


25 


805 


22.78 


1 .078 


i; 0 . 2 i; 


$0 


189 ii 


.88 


1 .091 


91 ;. 69 


% 


1996 


58.30 


1 .009 


99.79 



H. DISCUSSION. 

The sampling experiment suggests that the power of our sequential 
procedure is quite unsatisfactory against the particular Weibull alterna- 
tive examined for small N, but the power rapidly increases vjith N. Scru- 
tinizing of the density of (1|.2) indicates that there are relatively few 
"early failures", while a concentration of "v/earouts" probably accoiints 
for the rather long conditional expectation of the time to hit the 
boundary (reject the null hypothesis). 
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